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1 Introduction 



Stimulated by 't Hooft's H] and Polyakov's @ research van Nieuwenhuizen 
et al. || studied the influence of gravity on magnetic monopoles, but their 
results do neither guarantee the existence nor the stability of the solutions as 
shown in ||, ||, |J. Later Breitenlohner et al. @] and others (Lee et al. 0, 
Ortiz 0) investigated 't Hooft's and Polyakov's monopoles J2j) in curved 
space using numerical methods. Introducing the dimensionless parameters 
a = M w /gM Pl (with the Planck mass M P1 = 1/VG) and (3 = M H /M W , 
where My/ and Mh denote the mass of the corresponding Yang-Mills field 
and Higgs field respectively, Breitenlohner et al. Q presented the following 
results: For < (3 < oo the monopole mass varies only over a finite domain. 
Choosing the parameter /3 = solutions seem to exist for < a < a max . For 
a c ~ 1.386 < a < 1.403 ~ a max there are two solutions with different masses. 
Considering the mass as a function of a, one finds that it attains a maximum 
M maj! for a = a max , and that it decreases towards M P i/ 'g asymptotically as 
a tends to a c . This situation is shown in Figure 1. Each point represents a 
solution calculated numerically. The right plot of Figure 1 is a magnification 
of the left plot near the point B. The series of points marked by BC will be 
denoted by "upper branch" , the part AB by "lower branch" . 
In this paper stability of the nonabelian gravitating magnetic monopoles is 
investigated using the technique of spherically symmetric perturbations. Sta- 
bility considerations of similar nature have been conducted by Zhou and Strau- 
mann 0. The stability problem requires the numerical behaviour of the 
monopole solutions to be known. The construction of the monopole solu- 
tions in question (0]) is briefly reviewed. In order to investigate their stability 
the monopoles are subjected to small perturbations in time, which have a 
Fourier expansion due to its linearity. Eliminating the gravitational degrees of 
freedom, the remaining terms in the action are written as a quadratic form. 
Then the application of Jacobi's criterion (||, |0|) suggests instability on the 
branch BC and stability under spherically symmetric perturbations on the 
lower branch. Note however that this does not imply stability in general. In 
the final section the results are summarized and discussed. 



2 Solutions of the EYMH Equations 

The 't Hooft Polyakov monopole in curved space is a static finite energy solu- 
tion of an SO (3) gauge theory with a triplet of scalar fields. Using a coordinate 
free notation the action reads 

(1) 
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where R, F and $ denote the curvature scalar, the Yang-Mills field and the 
scalar field respectively. G and g are the gravitational constant and the gauge 
coupling. * is the Hodge star operator and Z7(3>) is a potential term which 
will be specified later. The gravitational field is described by a static spheri- 
cally symmetric metric tensor given by the line element ds 2 in Schwarzschild 
coordinates: 

ds 2 = A 2 (r)fi(r) dt 2 -^—dr 2 - r 2 (d9 2 + sin 2 9 d<j> 2 ). 

H(r) 

The minimal ansatz for the Yang-Mills potential A which is spherically sym- 
metric is denoted in polar coordinates and abelian gauge as follows (|Pfl,|jTU|, 

0): 

A = (W(r) T 2 + cot 9 T 3 ) sin 9 d<f> + W{r) T x d9. 

The Tj, i = 1, 2, 3, are the generators of the fundamental representation of 
SO (3). Using the potential A the field F is expressed by 

F = dA + ^ [A, A}. 

The Higgs field is chosen to be 

$ = H(r) (cos sin 9 T\ + sin sin 9 T 2 + cos 9 T 3 ) . 

In order to reduce the four dimensional action to a one dimensional integral 
the curvature scalar, the Yang-Mills, and the Higgs terms are substituted into 
the action ([I]): 

r°° / 1 

-tred 1 • • 1 



with 



S red = -j dr A[— (p + rp'-l)-^-^), (2) 



V 1 = W' 2 + jH' 2 , 

V 2 = ~ I)' + W 2 R 2 + lff2 r 2 (iy 2 _ 1)2- 

2r^ 8 

The term |/5 2 r 2 (if 2 — l) 2 represents the potential U (<§>). The equations of 
motion (EYMH equations) deduced from the action (0) contain the parame- 
ters a and (3 as well. Their values determine the stability properties of the 
monopoles. 

The "background solution" F(r) corresponding to global regular monopoles 
has been obtained solving the equations of motion with a multiple shooting 
method |]I2"]. 



2 



3 Stability Analysis 



The static background solution F(r) is subjected to a spherically symmetric 

perturbation f(r)e lujt , i.e. the perturbed and the background solutions are 
related as follows: 

F(r,t) = F(r) + f(r)e iuJt , (3) 

where F(r) is the vector of the functions H(r), W(r), Y(r) and M(r) and 
f(r) denotes the vector (h, v, y, m) of the radial parts of the perturbations. 
It is convenient to introduce M(r) and Y(r): 

Mr ) = 1 _^M ) and A{r) = e Y(r)_ 

r 

Instability in this context means the growth of the exponential term in time. 
The analysis is thus reduced to the question whether admissible perturba- 
tions, i.e. perturbations fulfilling the boundary conditions, with an imaginary 
frequency exist. Substituting F(r, t) into the action (HD one obtains a two 
dimensional reduced action differing from (Q) only by the time-dependent ad- 
ditional term 

v = -L (w 2 + -H 2 ) . 

An\ 2 ) 

The part of the reduced action S red quadratic in the fields and their derivatives 
turns out to be 

£red,2 = J°° dr Q Y f 2 _ ^ _ y \ _ ^ ^ + ^Jjf^ 
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(wl 2m 2W(W 2 - 1) 9 9 , 

+2y — + — Vx K — >-v - 2WH 2 v - 2W 2 Hh - —r H(H - l)h 

\a z r r z 2 

+ ^L{2W'v' + r 2 H'ti) - im {2v< 2 + r 2 h' 2 ) - 2 ^ W \ \~ X) v 2 - 2H 2 v 2 
- SWHvh - 2W 2 h 2 - ^-r 2 (3H 2 - l)h 2 - u 2 —{2v 2 + r 2 h 2 ) \ . 



Taking into account the boundary conditions at r = 0, the variation with 
respect to y yields 

fry fry 

-z + -zY' + D W ,L v' + D W L v + D H iL ti + D H L h = 0, (4) 

cr cr 

where L is the integrand of S red and Dw 1 , Dw, Dh>, and Dh denote the 
derivatives with respect to W, W, H' and H . (|) is solved by 

m = a 2 (2/iW'v + r 2 /j,H'h). 
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This result is used to eliminate m and y in S ' . Partial integration and 
introduction of h = rhj \/2 and the variable a by dr = Ajj, da leads to the 
normal form of a 2-channel Schrodinger problem S Ted ' 2 : 



5 rcd ' 2 = J™ da [D&DV + ^Q(r)^ - u 2 tfv} . 



(5) 



Here 



Q(r) 

Qu(r) 
Qia(r) 
Q 22 (r) 



/r \ nT f ,dh dv. 



Qn(r) Qi 2 (r) 
Q 2 i(r) Q 22 (r) 



^(h\V + a 2 \(A f irW'H' 
r V A 



Q 2 i(r) 



# 2 + 



- 1 



+ 4c^- 



1 AfxW' 2 , 



A K r 



The background solution is unstable, iff the Schrodinger problem admits a 
bound state, that is a solution with u 2 < 0. The standard technique for a 2- 
channel Schrodinger problem is to calculate the equations of motion from (|5|), 

Instead of 



which then define an eigenvalue problem with the eigenvalue 
analyzing the spectrum of this operator the positive defmiteness of ([5]) without 
the last term is investigated applying Jacobi's criterion ||, JT3|. Although in 
H it is formulated for compact intervals only the proof applies in this case 
to an infinite interval of integration, too, in view of the asymptotic behaviour 
of The discussion of Jacobi's criterion in || also provides a normalizable 
perturbation with negative energy. 

This criterion states that the monopole solutions are unstable iff 



det V(r) 



,(2) i 







for r G (0, oo). 



(6) 



hP~\ h.( 2 \ and denote the solutions of two linearly independent 

initial value problems for the equations of motion in the fields h and v 
(Jacobi equation). The initial values are chosen such that they are compatible 
with the symmetries of the background solution. 



A high order Runge-Kutta-Fehlberg method |12| has been employed for the 
numerical integration of this system of differential equations. 
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4 Results and Discussion 



First we discuss the results in the Prasad- Sommerfield limit. The numerical 
analysis shows that the branch BC in Figure 1 is unstable in the sense of 
linear perturbation theory, since the determinant has a zero in the interior of 
the integration interval. 

In Figure 2 the determinant (|6|) is shown for some values of the parameter 
a on the upper branch branch. For r Mm . n , that is the value of r where the 
function /i has a minimum (see H), detV(r) attains a local maximum. A 
second maximum at greater values of r is observed for some values of a. For 
a near a max the second maximum dominates, whereas, if a is approximately 
a c the left peak is well developed. For decreasing values of a (a — > a c ) the 
zero of the determinant moves to the left. 

It is convenient to rescale the determinant by fi(r)e~ ar in order to present 
the functions for all a -values of the upper branch in one plot (see Figure 3). 
The pronounced minimum in the left part of Figure 3 is due to the rescaling 
of detV(r). It suggests that lim Q ^ ac /x(r Mmin ) det ^(r Mmin ) = 0, i.e. in the 
neighbourhood of this point the determinant increases more slowly than 
For ct max the determinant has a zero at infinity in agreement with a theorem 
in |TJj which states that a change of the stability properties can only happen 
at a point of bifurcation in the masses. From this theorem it follows that the 
whole lower branch in Figure 1 consists of stable solutions, because of the 
stability of the Prasad- Sommerfield solution in flat space (a = 0). 
The bifurcation point is determined using Jacobi's criterion. For increasing 
values of j3 the method becomes more and more sensitive which means that 
the error becomes smaller and smaller. 

The determinant corresponding to solutions on the lower branch is plotted in 
Figure 4. As expected these functions have no zeros. The curve with the 
steepest ascent belongs to a = 0. 

Figure 5 and 6 show the results for some selected values of (3 ^ 0. One finds 
the same qualitative behaviour as for (5 — 0. In Figure 5 the determinant 
is plotted for given (5 and the corresponding maximal a-values. The lowest 
curve belongs to (3 = 0.05. Figure 6 uses the same /3-values, but a-values 
on the upper branch. Increasing f3 the zero of the determinant moves to the 
left. 

Some numerical results are collected in the following table. a c is obtained by 
linear extrapolation: 
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p 


Cmax 


^"max ( 


ru ■ ) 

Pram / 


M max 






0.00 


1.403034 


3.611 • 


10" 


-2 


1.000219 


1.3859 


0.05 


1.397555 


2.674 • 


10" 


-2 


1.000138 


1.3822 


0.10 


1.389459 


2.028 • 


10" 


-2 


1.000076 


1.3806 


0.15 


1.379974 


1.402 • 


10" 


-2 


1.000038 


1.3745 


0.20 


1.369676 


9.211 • 


10" 


-3 


1.000017 


1.3664 


0.25 


1.358883 


5.811 • 


10" 


-3 


1.000007 


1.3571 
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